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A /xm-scale calorimeter realized by using Schottky junctions as a thermometer is presented. Com- 
bined with a hybrid experimental method, it enables simultaneous time-resolved measurements of 
variations in both the energy and the heat capacity of subnanolitre samples. 



I. INTRODUCTION 

The present micromachining techniques allow the scal- 
ing of the calorimeter dimensions down to micrometer 
or even submicrometer scale, which results in high sen- 
sitivity and rapid response times. For instance, ther- 
mopiles fabricated on a thin SiN membrane have been 
used to determine the catalase activity within a single 
mouse hepatocytei Recently, nanowatt sensitivity and 
time constant of millisecond was achieved by optimiz- 
ing such a structured In general, membrane based mi- 
crocalorimeters using dc-methods have lately been un- 
der extensive study*^£&£ Additionally, ac-calorimetry 
methods 8 have been developed for measurement of heat 
capacity^ii^ and thermal conductivity^ of increasingly 
small samples. 

In this paper we describe the first /im-scale calorime- 
ter realized by using Schottky junctions as a thermome- 
ter. The ability to use present IC fabrication techniques 
makes Schottky junctions a very attractive choice for 
mass produced, low cost and high throughput calorime- 
ters for a variety of applications. Moreover, we utilize a 
novel measurement method which enables measurements 
of energy changes to be performed simultaneously with 
ac-calorimetry. 12 Combined with a short time constant 
(~20 ms) of the device, it allows direct time-resolved 
measurements of both the variations in the sample heat 
capacity and the energy changes due to for example a 
phase transition or other phenomena. The calorimeter 
may be used in ambient conditions and with liquid sam- 
ples allowing also real time monitoring, e.g., under a mi- 
croscope, and the operation is near isothermal utilizing 
only low heating rates, thus making it particularly suit- 
able for biological applications. The performance and 
reliability of the device and method were tested. We 
used subnanolitre drops of Dl-water as samples 2,4 and 
measured, simultaneously, the change of the heat capac- 
ity and the latent heat related to the evaporation of the 
drop. 



nm) SiN-layer (see Fig. [T]). The samples were fabricated 
on a lightly boron doped (12 Qcm) Si chip having a SiN 
layer on both sides. By using photolithography and Reac- 
tive Ion Etching (RIE), a square opening (1.2x1.2 mm 2 ) 
was etched on the SiN layer on one side. The opening 
was then used as a mask for chemical KOH etching to 
form a Si well. The etching was interrupted when there 
was still about 5-10 /im of Si at the bottom of the well, 
under which there was the SiN layer of the other side. 
After that the well was covered with SiN using PECVD 
deposition. This SiN layer was further e-beam patterned 
and RIE etched by using PMMA as a mask, to leave 
only a small SiN square at the center of the well bot- 
tom. A second KOH etching was applied to form a sep- 
arate Si island onto the SiN membrane. On the other 
side of that SiN membrane, openings on the SiN for the 
Schottky junctions to the Si island were done by e-beam 
lithography and RIE. The chip was exposed to chemical 
cleaning {procedure consisting of 2% HF (20s)/Piranha 
(5min)/2% HF (20s) }!£ prior to deposition of Ti, that 
acted as a metal for Schottky junctions. Ti was deposited 
by e-beam evaporator in UHV chamber. Ti layer was 
then patterned and (RIE) etched to form a heating ele- 
ment and the Schottky junctions onto the membrane. 




FIG. 1: A schematic and a micrograph (close-up) of the fab- 
ricated sample and the measurement principle. 



II. EXPERIMENTS 



A. Sample fabrication 

The fabricated calorimeter is composed of two Si/Ti 
Schottky junctions as a thermometer and a Ti heating 
element on a small Si island supported by a thin (600 



B. Experimental methods 

It is well known that Schottky diodes can be used as 
accurate thermometer at room temperature. The charge 
transport through Schottky junctions with low semicon- 
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ductor doping concentration is dominated by thermionic 
emission, which yields an exponential dependence be- 
tween the current and the temperature given by I t ot = 
AJst [exp(qV/k B T) - l]. 14 Here V is the voltage across 
the junction, A the area of the junction, q the electron 
charge, &b the Boltzmann constant and T the absolute 
temperature. The saturation current density Jst is de- 
fined by J ST = A*T 2 exp(-# jB //c B T), where A* is the 
effective Richardson constant and <\>b the barrier height. 
We consider temperature changes AT small enough so 
that the dependence of the current on temperature, J(T), 
can be linearized, i.e. J(T) ex (dlfdT)AT. In the present 
configuration there are two Schottky junctions in series, 
one being always forward biased and the other reverse bi- 
ased. The zero bias resistance of the fabricated samples 
with 35x45 /xm 2 Schottky contacts in ambient conditions 
is 5-10 and the variation of resistance about 10 %/K. 

At the equilibrium, the temperature of the Si island T 
is governed by the net power P in applied to it, and by 
the heat conductance K [W/K] of the SiN membrane to 
a heat bath formed by the Si chip at temperature To. In 
general, one can write 



T = T + P in /K, 



(1) 



if the variations in P^ n are slower than time constant of 
the calorimeter r = (C + AC)/K. Here C + AC is the 
total heat capacity of the Si island (C) and possible sam- 
ple on it (AC). An ac voltage V — Vo sm(uot) applied 
via the Ti- heater, i.e., the power Ph = Vq /Rx sin 2 (uot), 
where R is the resistance of the heating element, pro- 
duces by Joule heating a temperature fluctuation of the 
sample with a frequency 2uo. If uo ~ r _1 , the equilibrium 
result does not apply anymore and one obtains for the 
temperature of the island 



where the term AC describes the heat capacity of the 
sample which can also slowly vary in time, e.g., due to 
evaporation of a liquid sample drop. Thus, by simulta- 
neous measurement of the dc and ac signals, one is able 
to detect both the variation in temperature and the vari- 
ation in heat capacity, respectively. 

The simultaneous monitoring of the dc and ac signals 
was performed by using a measurement setup shown in 
Fig. [21 The dc-bias of the Schottky junctions was set to 
60 mV and the sinusoidal signal for the heating element 
had rms power of 37 fiW and a frequency of 10 Hz corre- 
sponding to the situation uo ~ r _1 . The dc component of 
the temperature sensitive current through the Schottky 
junctions was measured directly using DL Instruments 
low noise voltage preamplifier 1201 and to measure the 20 
Hz ac component, Stanford Research SR830 digital lock- 
in amplifier was used. These both signals were recorded 
by a computer with a sampling rate of 1 kHz. The fluctu- 
ating ac signal component was extracted from the mea- 
sured dc signal by digital notch filtering. 
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FIG. 2: A schematic of the measurement setup. 
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where up = arctan(2o;r) and Ps is the power produced 
by the sample (e.g. due to a chemical or biological re- 
action), or by some external source. Consequently, the 
dc-biased Schottky junctions will carry a current, directly 
proportional to T above, with one component having the 
frequency of 2uj and another being the dc component pro- 
portional to Ps an d the constant rms power fed in via the 
heater, VJ/2RK. 

Our calorimeter function is based upon the combina- 
tion of two calorimeter operation modes: the heat con- 
duction mode, which measures Ps/K via the dc compo- 
nent, and the ac-calorimeter mode with an ac voltage of 
frequency uo ~ r _1 applied to the heater (see Figs. [I] and 
[2]). According to Eq. (2), the amplitude of the produced 
2uo temperature component depends on the heat capacity 
of the sample as 



ST 



V 2 



1 



2R ^/K 2 + (2uo) 2 (C + AC) 2 



(3) 



C. Sample delivery 

The sample to be measured, i.e., drop of DI- water is 
delivered to the Si island of the calorimeter constructed 
onto the center of the SiN membrane (see Figs. [1] and [3]) 
using a very narrow pulled glass pipette. The injection 
is carried out by having the pipette back end attached to 
a pressure device (ASI MPPI-2 microinjector) which can 
supply a short pressure pulse with adjustable pressure 
and pulse time. In principle, this device allows one to 
precisely and reproducibly (from measurement to mea- 
surement) control the amount of liquid injected. How- 
ever, it is not possible to quantify the amount since the 
total amount depends on the viscosity of the liquid, the 
surface effects between the liquid and the pipette, the 
pipette dimensions, etc. 

All the time during the sample delivery and the mea- 
surement, the end of the pipette is at the very near prox- 
imity of the Si island of the calorimeter. This is necessary 
since, due to high surface tension of the water, the small 
drops used in this experiment will not have sufficient 



mass to detach from the tip of the pipette unless they 
touch the Si island surface. During the pressure pulse, 
the drop size increases until it touches the hydrophilic 
surface of the Si island and spreads on. After the pulse, 
the drop size starts to decrease due to evaporation of the 
drop and also due to the reflux of the water back to the 
injector as will be discussed later. 

After the evaporation of the drop, a small amount of 
water is left to thermally bridge the pipette and the Si 
island (see Fig. [3|). This is due to hydrophilicity of the 
surface of the Si island and the surface tension of wa- 
ter, and it creates an additional pathway for heat trans- 
fer between the calorimeter and the environment. Since 
the relative positioning of the tip of the pipette and the 
calorimeter remains constant throughout the whole ex- 
periment, it is expected that this additional pathway for 
heat transfer will simply increase the total heat conduc- 
tivity K in our simple model by a constant factor which 
can be experimentally determined. Thus the total ther- 
mal conductivity may be assumed to be constant over 
the time span of the experiment. 

Prior to measurement of the device parameters (see 
next section), a drop of water was injected onto the Si is- 
land surface which was then allowed to evaporate to form 
the water bridge between the calorimeter and the pipette 
as described above. This configuration is now considered 
to be the empty device, since after each injection the sys- 
tem returns to this state. The pipette is kept on its place 
during all the further experiments. Thus, the obtained 
parameter values already include the effect of the addi- 
tional heat conductivity through the water, i.e., the effect 
of the pipette slightly touching the device. Importantly, 
this effect can be accurately measured, and it remains 
constant through the whole measurement process, which 
makes it controllable. However, for more practical use 
of the device, other delivery systems, e.g., microfluidistic 
channels, could be implemented. 



D. Calibration 

The dc signal (with 60 mV dc-bias) dependence on 
the temperature was calibrated by using a Pt-100 ther- 
mometer as a reference, resulting in a linear relation with 
a slope of 9.5 x 10 -10 A/K between the current and the 
temperature near the operation temperature 23° C. Ad- 
ditionally, the relaxation method was used to determine 
the thermal conductance, and heat capacity, C, of the 
empty device (including the pipette for water, see above 
and Fig. [3|) without the sample: To obtain the thermal 
conductance according to the Eq. (pQ), we measured the 
dc current as a function of slowly increasing dc power (ap- 
plied to the heating element), yielding 6.95 x 10 -5 W/K. 
The time constant of the calorimeter was determined by 
measuring the dc signal response to the step pulse fed 
into the heater, resulting in a time constant of 21 ms. 
Finally, the calorimeter heat capacity was obtained as 
the product of the time constant and the conductance 
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FIG. 3: (a) - (d) Cartoon of evaporation of a water drop after 
injection. In the end, a small amount of water stays to bridge 
the end of the pipette and the Si island of the calorimeter. 
This configuration, shown in (d), is considered as the "empty 
device" . The thermal conductivity K and the heat capacity C 
are determined for this "empty device", thus accurately and 
systematically including the effect of the pipette. 

which yielded a value of 1.53 x 10 -6 J/K. 

III. RESULTS AND DISCUSSION 

A. Measured signals for subnanolitre drops of 
water 

The validity of our measurement method was evaluated 
by injecting several subnanolitre-scale drops of DI- water 
with different (undetermined) volumes onto the calorime- 
ter, and by measuring the ac and dc responses as func- 
tions of time as explained in section III Bl In Fig. [4^) is 
shown the dc and ac signals of the thermometer with sev- 
eral injection times of the nanoinjector, ranging from 10 
ms up to 60 ms (injection pressure being constant ~ 0.5 
Bar). During all the measurements, the pipette was in 
thermal contact with the device as described above in 
section III CI 

The figure 0t>) illustrates the function of the sensor in 
detail for 60 ms injection time. From the figure one can 
see the abrupt change in the two signals when the injec- 
tion of drop takes place, i.e., at around 288.5 s. The dc 
signal increases, indicating the temperature decrease due 
to evaporation of the drop while the ac signal decreases 
due to increase in heat capacity of the system. Between 
the time 289-295.5 s the dc signal remains approximately 
constant, which implies that the power consumption due 
to evaporation (Pg in Eq. fl2J)) and thus the surface area 
of the drop, remains almost constant. This is due to hy- 
drophilicity of the Si island surface, which causes a flat 
shape of the drop and the height of it to vanish almost 
completely before the horizontal area starts to decrease 
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during evaporation (see Fig. [3]). This also implies that the 
volume of the evaporated part of the drop is decreasing 
with approximately constant rate. During this timeframe 
the ac signal is increasing which is due to continuous de- 
crease in heat capacity, AC, of the system as the drop 
evaporates, which is also shown as a dash-dotted (green) 
curve in Fig. SJd). As the heat capacity is proportional 
to volume, this is consistent with the conclusions made 
from the dc signal. 

In figure UJd) the dashed line shows the calculated 
(along Eq. (3)) variation of ac signal amplitude as a func- 
tion of time assuming a constant reducing rate of volume 
after the injection. Also, a finite time for injection, dur- 
ing which the drop is growing linearly, is assumed in cal- 
culation. These assumptions can be clearly verified from 
the curve showing AC as a function of time. Heat ca- 
pacity of the drop, AC, as well as the heat capacity, C, 
and conductance, K, of the calorimeter for the calcula- 
tion are obtained from the measurement and calibration 
data. The excellent agreement of the calculated curve 
with the measured one further verifies that the reducing 
rate of the total heat capacity of the sample during the 
evaporation is approximately constant in time, as the dc 
signal already suggested for the evaporated part. Since 
the total reduction of the heat capacity is due to evap- 
oration and also due to the reflux of the injected water 
back to the pipette, as will be further explained in next 
section, this implies that also the reflux happens at a 
constant rate. 
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FIG. 4: (a) The dc (blue) and ac (red) signals of the ther- 
mometer for several consequent injections, the injection times 
ranging from 10 ms up to 60 ms. (b) A closeup of both sig- 
nals with 60 ms injection time. The dash-dotted (green) line 
is the heat capacity of the drop, AC, obtained from the ac 
signal as a function of time and the dashed line (black) shows 
the theoretical signal calculated along Eq. (3) assuming linear 
decrease in AC. The two arrows show the corresponding axes 
for the curves. 



B. Detection of nonidealities in sample delivery 

In general, the heat capacity, AC of the sample can 
be expressed in terms of specific heat c and mass m of 
the sample as AC = cm. Furthermore, for samples that 
evaporate, AC = cm = cE/£, where E is the total energy 
required for evaporation of a drop and £ is the latent 
heat of evaporation. From the measured dc signal one 
is able to obtain the energy E as the time integral of 
the signal, and the variation of the heat capacity of the 
system is related to the ac signal via Eq. ([3|). Such com- 
plementary information allows to understand details of 
the process and eliminate the effect of measurement non- 
idealities (which are of increasing importance when the 
sample volume decreases), as will be shown in the fol- 
lowing. Furthermore, for evaporating samples, it offers a 
test of consistency since the heat capacity can be deter- 
mined in two ways: from the ac signal, and, in addition, 
from the dc signal via the relation AC = cE/i. 

As the main non-ideality in the measurement, it was 
observed that the non-ideal backpressure applied to the 
nanoinjector produced a reflux of the injected drop back 
into the pipette during evaporation (the backpressure is 
used to compensate the capillary forces that result in a 
reflux of the injected liquid back to the pipette). There- 
fore the energy of evaporation E obtained by integrat- 
ing the dc signal does not give, without a correction, 
the initial heat capacity of the sample from AC = cE/£. 
The required correction was obtained by performing mea- 
surements for drops of different (unknown) volumes until 
they were fully evaporated, and by 1) determining the 
mass of the initial drop from the ac signal as the ratio 
mc = AC/c of the maximum variation of heat capac- 
ity AC and the specific heat of water c, 2) determining 
the mass of the evaporated part of drop from the dc sig- 
nal as the ratio tue = E/£ of the energy required for 
evaporation E and the latent heat of water £. The dif- 
ference between the mass of the initial drop and that of 
the evaporated part then gives the mass of the part that 
was refluxed to the injector. 

A plotting of mc against tue for each drop, as shown 
in Fig. O reveals that the mass tue obtained from evapo- 
ration energy is systematically smaller than the mass mc 
obtained from minimum of the ac signal. This result is 
easy to understand by noting that since the pipette sucks 
some of the liquid back during the evaporation, the mass 
ttie obtained from evaporation energy must be smaller 
than the mass mc obtained from the initial heat capac- 
ity of the drop (ac signal). Furthermore, since the reflux 
is due to a constant pressure difference, it is supposed to 
happen at approximately constant rate, which can also 
be verified from the signals in Fig. 0] as explained in pre- 
vious section IIIIBI As shown by the Fig. as much 
as 72% of the initial mass is lost by the reflux in every 
pipetting. 

Such measurements can also be used to calibrate 
the pipette for different liquids, surfaces, etc., which 
is often found to be a problem in nano- and subnano- 
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litre injections. Therefore, the hybrid method and 
the device we use allow to accurately distinguish vari- 
ous non-idealities of the measurement from the chem- 
ical/physical/biological phenomena of interest. This is 
certainly a strength, especially at nanoscale where non- 
idealities in the measurements are difficult to avoid. 
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FIG. 5: Relation between the initial mass of the drop mc and 
the mass of the evaporated part the- Circles correspond to 
the measured data and line is a linear fit to it. 



C. Evaluation of the calorimeter 

Finally, as a test of performance and consistency of our 
method, we present in figure [6] the following analysis: We 
performed measurements for different sizes of the drops. 
For each drop, the information given by the two signals, 
ac and dc, is plotted as a point in a two dimensional graph 
where one axis is the minimum of the normalized ac lock- 
in signal (corresponds to the heat capacity of the initial 
sample via Eq. (j3j)), and the other axis is the integrated 
dc signal giving the evaporation energy E. The latter 
yields the mass of the evaporated part of the drop the 
via the relation = Ejt. However, we cannot use the 
mass mE as such since it excludes the part gone back to 
the pipette by the reflux, and when considering e.g. the 
minimum value of the ac signal, one has to use the mass 
of the initial drop, which includes the part that later will 
be gone back to the pipette by the reflux. From the 
analysis presented above, we know the ratio between the 
evaporated and the initial masses to be constant, ~ 0.28. 
Thereby we can use a corrected value of mass, i.e. m = 
mc ~ raE/0.28 = E/(0.2S£) and thus corrected energy 
Eq = .E/0.28, which is plotted as x-axis in[6l 

As a fit to the theory and test of consistency, we 
then take the measured values of E, use the relation 
AC = cm = (c/£)Eq (with c/i as the only fitting pa- 
rameter), insert this value of AC in Eq. ([3]), and plot the 
resulting curve in the figure. It corresponds excellently 
to the measured data points. The values for C and K 



used in this analysis are taken from the measured cali- 
bration data. From the value given for c/i by the fitting 
and using the known value for c or i, we obtain the cor- 
rect value for i or c, respectively, with the accuracy of 
about 7%. This accuracy characterizes the overall perfor- 
mance of the device and method at this stage, and can be 
improved by optimizing the sample fabrication and the 
measurement setup. Note that for evaporating samples 
where one of the parameters c, i or volume is known, the 
method can used for determining the two others. 

The ac signal also allows high accuracy mass measure- 
ments: the smallest measured change of heat capacity 
(for 10 ms injection time) was approximately 10 _T J/K 
which, for water, corresponds to 20 ng mass. This could 
be extremely useful for example when measuring a re- 
action heat of some chemical or biological reaction, since 
one is able to simultaneously measure the amounts of the 
reagents added and the energy produced or consumed by 
the reaction. 
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FIG. 6: The measured minimum values of the normalized ac 
signal and the measured total energy required for drop evap- 
oration (with the reflux correction), together with the theo- 
retical fit obtained using the relation of evaporation energy 
and heat capacity and Eq.©. 



IV. CONCLUSIONS 

In summary, we have developed a microcalorimeter 
capable of providing simultaneous quantitative time- 
resolved information of the sample temperature, power 
and total energy production/consumption as well as of 
sample heat capacity and its variations. The small heat 
capacity of the device which enables fast and sensitive 
measurements carried out in room temperature, together 
with the mass production potential of the Schottky junc- 
tions, hold a promise for high throughput and cost effec- 
tive tool for biological and chemical applications. 
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